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Abstract 

The observed acceleration of the universe, explained through dark energy, could alternatively be explained through a 
modification of gravity that would also induce modifications in the evolution of cosmological perturbations. 
We use new weak lensing data from the COSMOS survey to test for deviations from General Relativity. The depar- 
ture from GR is parametrized in a model-independent way that consistently parametrizes the two-point cosmic shear 
amplitude and growth. 

Using CMB priors, we perform a likelihood analysis. We find constraints on the amplitude of the signal that do not 
indicate a deviation from General Relativity. 

Key words. Gravitational lensing. Cosmology: Large-scale structure of Universe, Cosmological parameters. 



1. Introduction 

The ACDM model is currently the best fit to the available 
cosmological data. In this framework, the period of accel- 
erated expansion of the Universe is due to dark energy, a 
dominant component of the energy density with a negative 
equation of state. A possible alternative, worth to inves- 
tigate observationally, is to replace the standard gravita- 
tional theory of General Relativity (GR) by a modification 
that admits self-accelerating solutions. Such gravity theory 
must converge to GR both at high redshift, when the ex- 
pansion is not accelerated, and on small scales, since GR 
has consistently passed local tests of gravity (jWill 2006|) . 

The various classes of modified gravity and dark en- 
ergy models (jUzan 2007[) cannot be distinguished based 
only on observations of the cosmological background. They 
are only distinguishable due to their different predictions 
for cosmological structure formation ( |Zhang et al. 2007] ). 
The evolution of cosmological background and perturba- 
tions have been studied in various classes of physically mo- 
tivated modified gravity theories. Most prominent today 
are theories where the modifications arise from extra di- 
mensions, with matter fields confined to a 4-dimensional 
brane embedded in a higher dimensional bulk, such as the 
Dvali-Gabadadze-Porrati (PGP) mo del (jDvali et al. 20001 
|Deffayet et al. 200 2 Luc et al. 2004 ), and also the s o-called 
i(R) theories (jCarroll et al. 20041 |Song et al. 2007] . In the 
latter, the field equations are modified because the corre- 
sponding Einstein action no longer depends linearly on the 
Ricci scalar (R) but there are extra terms on R. 
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The standard way of studying cosmological struc- 
ture formation is to introduce scalar perturbations 
in the Friedmann-Robertson- Walker (FRW) space- 
time metric. In the Newtonian or longitudinal gauge 
(Mukhanov et al. 1992), where by construction the metric 
remains diagonal, there are only two scalar perturbations: 
the scalar fields ^ (the Newtonian potential) and $ (the 
spatial curvature potential). They are sourced by the cos- 
mological fluid perturbations contained in the stress-energy 
tensor. Conservation of the stress-energy tensor imposes 
evolution equations for its components. In addition, the 
equations of the theory of gravitation, relating the metric 
components with the energy components, bring extra 
evolution and constraint equations. 

In modified gravity, structure formation has been stud- 
ied with various parametric and non-parametric approaches 
( |Bertschinger fc Zukin 20081 ISilvestri fe Trodden 2009)) . 
Parametric approaches include direct parametrization of 
the linear growth (jLinder 20051 ILinder fc Cahn 2007j) : 
the so-called parametrized post-Friedmann frame- 
work (PPF) (|Hu fc Sawicki 2007|) and several related 
phenomenological post-GR parametrizations, which 
give a model-independent description of deviations 
to GR in a way that is convenient for the test- 
ing of cosmological predictions (jCaldwcll et al. 2007; 

lAmendola et al. 20081 IDaniel et al. 20081 

IDaniel et al. 2 010; Pogosia n et al. 2010[ ); and phenomeno- 
logical applications to parametrizations of particular 
theories of modified gravity ( |Acquaviva et al. 2008; 
IZhao et al. 20091 |Song et al. 2010[ ). 

The various metric and energy perturbation variables 
in modified gravity scenarios may be measured by differ- 
ent probes (|Huterer fc Linder 20071 |Jain fc Zhang 2008 
ISchmidt 200j |Song fc Pore 20091 IGuzik et al. 2010 
Bcy non et al. 2010[ ). In particular, weak lensing by large- 
scale structure (cosmic shear), which probes the sum <fr + \E r , 
is well suited for such studies (jUzan fc Bcrnardcau 2001; 
IKnox et al.~20 06). In cosmic shear measurements, the 
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gravitational information is contained in the equations 
relating the measured matter density perturbations with 
the two potentials, or equivalently in one of those equa- 
tions (a generalized Poisson equation) and an equation 
relating the two potentials (the anisotropy equation). A 
parametrization of those two relations is thus well suited 
for cosmic shear measurements, which will provide a 
degenerate estimation of those two quantities. This is the 
approach we follow in this work. 

We compute cosmic shear two-point correlation mod- 
els in tomographic bins, using standard cosmological 
parameters plus two modified gravity parameters, and 
analyse them against the new COSMOS tomographic 
cosmic shear measurements of Schrabback et al. (2010). 
Other cosmic shear data have already been used in recent 
studies of modified gravity, using similar or different 
parametrizations (jDore et al. 20071 IThomas et al. 2009 
IDaniel et al. 20091 



Zha o et al. 20T0l 



Bean fc Tangmatitham 2010[ 
IDaniel fc Linder 20101) ! 



The stan- 
dard parameters are constrained by using WMAP 7-year 
CMB data (jKomatsu et al. 2010|) . Our goal is to study the 
impact of this phenomcnological model in the degeneracy 
between the amplitude of the cosmic shear signal and the 
linear growth of structure. In other words, to investigate 
if, given the strong cosmological priors imposed by the 
WMAP7 results, there is room to compensate a possible 
alteration of the predicted redshift-dependent amplitude of 
the cosmic shear effect with a modification of the growth 
of matter perturbations, which, phenomenologically, would 
keep open the possibility of an alternate gravity theory as 
an origin for an effective dark energy behaviour. 

In the next section, we present the parametrization of 
modified gravity used. The new COSMOS dataset is de- 
scribed in Sect. 3. The statistical analysis and its results 
are presented and discussed in Sects. 4 and 5. We conclude 
in Sect. 6 with a reminder of the main assumptions made. 



2. Parametrization 

Images of background galaxies are distorted (by conver- 
gence and shear) in the perturbed FRW spacetime. In the 
Newtonian gauge, considering scalar perturbations and a 
flat background, this spacetime is described by 

ds 2 = - (1 + 2V(x, t)) At 2 + a 2 (t) (1 - 2<f>(x, t)) S ij dx i dx j . 

(1) 

The local deflection of null geodesies, with respect to un- 
perturbed ones, depends on the traveling time of a light 
ray, and thus on the sum of the metric perturbations, 



^ « 1 -(* + $) . 
dx 



(2) 



Integrating the multiple deflections over cosmological dis- 
tances, and introducing the lens equation, the cosmic con- 
vergence is written, to first order in the amplification matrix 
and for a flat Universe, as (Bartchnann & Schneider 2001) 



k(6, w) 



2 Jo 



, . w — w ,~ 



w'V 2 + $] (w'O, to') , (3) 



where (9, w) are comoving spatial coordinates and the 
Laplacian of the lenses potentials at comoving distances 
w' is integrated up to the source at distance w. 



The equations of the theory of gravitation at perturba- 
tion level relate the metric perturbations (the potentials , 
$) to the perturbations of the cosmological fluid (matter 
density contrast S, pressure perturbation 5p, fluid veloc- 
ity v, and stress or anisotropic pressure a). In GR these 
equations are the linearized Einstein equations. The phe- 
nomcnological approach parametrizes deviations from GR 
by introducing extra functions in the linearized Einstein 
equations. 

Following Ma & Bertschinger (1995), a combination of 
the — and — i linearized Einstein equations relates $ 
with the comoving matter density perturbation: it is the 
generalized Poisson equation. In Fourier space, with k 2 be- 
ing the spatial Laplacian, it writes, 



fc 2 $ = -4wGa 2 pA 



317, 



>H 2 A 



(4) 



where A is the comoving density perturbation A = 5 + 
3H(1 + w)v/k 2 , w is the equation-of-state parameter w — 
p/p, and the second equality is written for dark matter 
inserting f2 m . The deviation from GR is parametrized by 
introducing the mass screening phenomenological post-GR 
function Q(a, k), replacing the gravitational constant by an 
effective function, 



G cff = GQ => fc 2 $ = -A-k G Q a 2 pA . 



(5) 



Notice that a hypothetical clustering dark energy compo- 
nent would also add a contribution to the Poisson equation. 
Factorizing out the dark matter density contrast, a multi- 
plicative term, function of <5de/<5> would appear, and play 
the same role as Q. 

The i — j Einstein equation gives a second constraint, 
relating the difference <f> — W with the stress a: it is the 
anisotropy equation, 



k 2 (<t> - *) = 12tt G a 2 p (1 + w)a. 



(6) 



In GR the two potentials arc identical. The deviation from 
GR is parametrized by introducing the gravitational slip 
phenomenological post-GR function rj(a,k), 



* = (1 + 77)$ 



/c 2 ($ 



AirGQr 1 a 2 pA. (7) 



Again, a hypothetical dark energy component with 
anisotropic stress a would have the same effect as Qr\. 

The two remaining independent Einstein equations, for 
example — i and i — i, are evolution equations for and 

and will not be needed when considering the quasi-static 
approximation, valid for sub-horizon scales. 

The linearized Einstein equations allow thus for 2 inde- 
pendent post-GR functions. From Eqs. (J5J) and (0, we can 
also write modified equations for the Laplacians of $ + ^ 
and They are 



/c 2 ($ + *) = -SttGQ (1 



a 2 pA 



(8) 



and 



fc 2 * = -4tt G Q(l + 77) a 2 pA . 



(9) 

These equations introduce two derived post-GR functions 
that are useful for cosmic shear studies: £(a) = Q(a)(l + 
ry(a)/2) and T(a) = Q(a)(l + 77(a)). 

One approach, taken e.g by Daniel & Linder (2010) 
and Zhao et al. (2010), to study the post-GR modifica- 
tions as function of scale and redshift is to parametrize 
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them with various constant parameters in redshift and scale 
bins. Here, we consider them to be scale-independent in the 
limited range of scales provided by the COSMOS data, al- 
though it is clear that in a broad range of scales there is 
scale dependence at least to differentiate the GR behaviour 
on solar system scales from the possible modifications on 
cosmological scales. 

Regarding the redshift dependence, we adopt a func- 
tional form using the scaling argument of Caldwell et al. 
(2007), according to which the departure from GR orig- 
inates an effective smooth dark energy A with constant 
w = —1. The deviations in the growth history scale with 
the deviations in the expansion history, such that they do 
not turn on too early, deviating from GR at early times, 
or too late. In this approach, the background evolution is 
identical to ACDM, no further modification is assumed in 
the Friedmann equations, and the post-GR functions take 
the GR values at a < 1 (Q = 1, r) = or E = 1, T = 1) 
and deviate from GR proportionally to pde{o-) / Pm(a) ■ This 
implies that 



Q{a) = l + < 



77(a) = 7i a a 3 



(10) 



This dependence for the fundamental functions Q(a) and 
77(a) implies a different redshift dependence for the derived 
functions £(a) and T(a). Only in the case of small parame- 
ter values, (Qa < 1, ?7a < 1), do the derived functions have 
a similar scaling, 



£(a) 



1 



I» 



(11) 



where £ a = Q & + 77 a /2 and r a — Q a + ?/ a . We will work 
in this ansatz, which reduces the range of values allowed 
for the post-GR parameters. In particular, we will consider 
models such that the difference between the two potentials 
is less than 50% and the effective gravitational constant 
does not differ more than 50% from the GR value. Large 
deviations are better studied with direct parametrizations 
in redshift bins, in the absence of a well-motivated func- 
tional form. Also, the approach of modifying the linearized 
Einstein equations should break down for too large devia- 
tions. 

Going back to the convergence, Eq. (|3|), we co mpute its 
power spectrum (Bartelmann & Schneider 2001) and, in- 
serting Eqs. ([5]) and ([7]), obtain the cosmic shear power 
spectrum in the following form, 



a 2 (w) 



(12) 

The amplitude of the power spectrum depends on the de- 
rived function E(a). This equation considers the correlation 
of source galaxies from two redshift bins, with redshift 
distributions p (w) , which are lensed with geometrical effi- 
ciency factors given by, 



9i(w) = / dw'pi(w') 



s w' — w 
w' 



(13) 



Deviations to GR thus modify the cosmic shear ampli- 
tude via S. They further impact the cosmic shear power 
spectrum via a modification of the matter power spectrum, 
Ps, present in Eq. (IT21) . The conservation of stress-energy 
provides the evolution equations for the matter perturba- 
tions. The evolution equations for S and v depend on the 



stress-energy components, the time evolution of the po- 
tentials and on the Laplacian of the Newtonian potential 
(Ma & Bertschinger 1995). Inserting the constraint equa- 
tions Eqs. (J5J and ([7]), and writing for the comoving den- 
sity A, a great number of terms arise. The solution is a 
scale-dependent growth in the linear regime, even for scale- 
independent post-GR parameters (Daniel et al. 2010|) . 

For the case of sub-horizon scales, assuming time deriva- 
tives of the perturbations are negligible compared to the 
spatial derivatives, and A « 8, the evolution equation of 
the cold dark matter perturbation greatly simplifies. It only 
contains a Hubble drag term and a term on k 2 ^> which 
means, from Eq. that the growth depends directly on 
the function T(a). 

In this regime, the growth solution can be cast in a 
form that explicitly separates the expansion history and 
gravitational dependencies of the growth (jLinder 2 005). 



dln(<S/a) 
dln(a) 



(14) 



This introduces the growth parameter, which contains the 
gravitat ional dependence, a nd has the value 7 = 6/11 
in GR ( |Wang fc Steinhardtp . In the case of ACDM back- 
groun d evolution and scaling evolution of T(a), it takes the 
value (|Linder fc Cahn 20071 IDaniel et al. 2010[) , 



7 = TT 



1 - 



2 1 - fL 



(15) 



3. COSMOS cosmic shear data 



The new COSMO S cosmic shear catalogue 
(Schrabback e t al. 2010}) (hereafter S10) contains 
around 450 000 galaxies with magnitude i < 26.7 in 
an area of 1.64 deg 2 , corresponding to a density of 76 
galaxies /arcmin 2 . 

The cosmic shear analysis of S10 used the CO SMOS-30 
photometric redshift catalogue (|Ilbert et al. 2009)) . which 
provides individual photometric redshifts for galaxies of 
magnitude i < 25. Those were grouped in 5 redshift bins. 
The highest magnitude galaxies (i > 24) of the lowest red- 
shift bin (z < 0.6) suffer from catastrophic errors due to 
the dificulty to locate the 4000 A break in their spectra. 
They are thus likely to contain high-redshift galaxies and 
were removed from the analysis. In addition, the redshift 
distribution of galaxies with i > 25, for which individual 
redshifts were not available, was estimated in S10 and those 
galaxies were placed in a single wide redshift bin. Table [T] 
shows the 6 redshift bins used, containing the remaining 
roughly 420 000 galaxies. 

The data used in this work consists on a vector of shear 
correlations and £_ computed with the estimator given 
in Eq. (13) of S10. We use all the cross-correlations between 
redshift bins (12, 13, 14, 15, 23, 24, 25, 34, 35, 45, 16, 26, 36, 
46, 56) and the auto-correlation of the wide sixth bin. The 
auto-correlations of the narrow bins (11, 22, 33, 44, 55) were 
not used to minimize the effect of intrinsic alignments. In all 
cases, the shear correlations were computed on 5 angular 
separation bins ranging from 0.7 to 21 arcmin. The data 
vector has thus (15 + 1) x 5 x 2 = 160 components. S10 
found the B-modes to be consistent with zero at all scales. 
In this case, we work directly with the £± correlations and 
there is no need to make an E/B mode decomposition of the 
shear field, which has the disadvantage of being non-local. 
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Table 1. Description of the redshift bins used in the anal- 
ysis, with their amount of contributing galaxies and mean 
redshifts. From Schrabback et al. (2010). 



Bin 


-2 m in 


-2 m ax 


TV 


<^> 


1 


0.0 


0.6 


22 294 


0.37 


2 


0.6 


1.0 


58194 


0.80 


3 


1.0 


1.3 


36 382 


1.16 


4 


1.3 


2.0 


25 928 


1.60 


5 


2.0 


4.0 


21718 


2.61 


6 


0.0 


5.0 


251958 


1.54 ±0.15 



0.0 t 

0.6 



1.0 

Q 



0.0 t 

0.6 



1.0 
I+77 



Figure 2. Marginalised one-dimensional PDFs obtained for 
the phenomenological post-GR parameters. 
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Figure 1. Correlation function £S_ 6 . Top panel shows the 
measurement (points with error bars), the best- fit GR 
model for a non-tomographic analysis (O m = 0.27, <r s = 
0.76) (solid), modified gravity models with the same stan- 
dard parameters but (E = 1.6, 7 = 0.55) (dashed), or 
(E = 1, 7 = 0.65) (dotted). Bottom panel shows the de- 
viation between the same dashed and dotted models with 
respect to the GR model. 



The 160 x 160 covariance matrix of the shear correlations 
was estimated in S10 from 288 realisations of COSMOS- 
likc fields obtained from ray-tracing through the Millenium 
simulation (jHilbert et al. 2009|) . The inverse covariance ma- 
trix, needed for the likelihood analysis, was debiased by 
applying the correction of Hartlap et al. (2007). 

4. Analysis 

We build models on a three-dimensional grid consisting of 
(Q a ,rj a ) and a redshift-calibration parameter f z . 

The parameter f z is a nuisance parameter that allows to 
include an uncertainty on the redshift calibration of bin 6. 
Following S10, we marginalize over this uncertainty by in- 
cluding shifted redshift distributions on the grid of models, 
such that each model is evaluated assuming various redshift 
distributions p (f z z) in the range 0.9 < f z < 1.1. 

We choose to work with the two original phenomenolog- 
ical parameters, to be able to have an uniform coverage of 
the (Qa,77a) plane. We probe a relatively narrow range of 
values: —0.5 < Q a < 0.5 and —0.5 < ?/ a < 0.5. As discussed 
in Sect. 2, this choice allows the redshift-dependence of the 
post-GR derived parameters to scale with the background 
expansion. The grid mapped on the plane of the derived 
parameters, the ones that have a direct impact on the cos- 
mic shear signal, i.e., the amplitude E a and the growth T a 
or 7, introduces a correlation between these 2 parameters. 



Its range is such that E(a = 1) can deviate from the GR 
value, E(a = 1) = 1, by up to 70%, while 7 extends to a 
maximum deviation of 40% of the GR value of 7 = 0.55. 
(Hereafter we denote the values of the post-GR parameters 
at a = 1 by Q, r), E and T). 

The grid is computed for each point of a sample of stan- 
dard cosmological models built by thinning-out the ACDM 
WMAP-7 year Monte Carlo Markov Chains available in the 
WMAP archive B In this way, the analysis has strong pri- 
ors on the standard cosmological parameters. The average 
as and ft m values of the models used are as — 0.807 and 
n m = 0.271. 

For each model, we evaluate the various shear power 
spectra using Eq. (fT2"j) . The matter power spectrum is com- 
puted using the transfer function of Eisenstein & Hu (1998) 
and the linear growth of Eq. flT4"]) with 7 given in Eq. (fT5|) . 
Since the small scales of the COSMOS data are in the quasi- 
linear and non-linear regime, we need a non-linear prescrip- 
tion to correct the linear power spectrum. We compute the 
non-linear power spectrum using the formula of Smith et al. 
(2003), which was built from N-body simulations of ACDM. 
Finally, the shear correlations are computed from 

CI/S0) = ^l°° d££J 0/4 (te)P^(£) , (16) 

where J„ is the n th order Bessel function of the first kind. 

As an example, Fig. [T] shows the auto-correlation 
function for the wide redshift bin. In the models shown, 
E and 7 take the extreme values used in the grid, when 
the other parameter is kept fixed at GR. Consistently with 
the approximations described before, the effect is scale- 
independent. We notice the impact of 7 is very small on the 
COSMOS angular scales whereas the impact of E clearly 
dominates. The effect of the growth parameter would be 
more important on larger scales, due to the k 2 factor in the 
Poisson equation (see also lDaniel et al.~2 010). 

The likelihood of each model of the grid is evaluated 
with respect to the previously described COSMOS data 
vector of dimension 160. The likelihood values are then used 
to update the original weights of the chain, in order to 
obtain the joint CMB-cosmic shear weights according to 
the method of importance sampling (Lewis fc Bridle 2 002). 

The resulting marginalised Probability Distribution 
Functions are shown in Fig. [2j where it is clear that the 
phenomenological parameters Q and rj are not constrained 
by the data in the range of values probed. 
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Figure 3. The dotted lines show the limits of the probed 
region. The 68% C.L. likelihood contour essentially fills up 
that region, showing that for these data the range of models 
chosen constitutes a significant prior. 
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Figure 5. Marginalised likelihood contours on the (E,<7g) 
plane at 68% and 95% confidence levels, obtained from the 
full analysis (filled contours), and using only the '66' cor- 
relation (open, dashed lines). 
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Figure 4. Marginalised one-dimensional PDF of £. 
5. Results and Discussion 

We consider now the derived parameters, E and 7, the ones 
that directly affect the cosmic shear signal. Fig.|3]shows the 
result in the (2,7) plane. In this analysis, given the scale- 
independent approximations made (both in the choice of 
scale-independent parameters and in the use of small scales 
where the modifications to growth are scale- independent), 
the post-GR models differ only in their amplitudes and 
redshift-evolution of the amplitudes. E and 7 are thus de- 
generate and we constrain the following combination of pa- 
rameters: 

E + 2.3 (7 - 0.55) = 0.99 ± 0.31 (68%C.L.). (17) 

This result combines information from the data, the cos- 
mological prior, and the prior given by the range of the grid. 
The latter introduces a correlation between £ and 7 that 
arises from the fact that the uncorrelated grid on the (Q, -q) 
plane of the phenomenological parameters maps onto a cor- 
related grid on the (E, 7) plane, as shown in Fig. [3J 

The corresponding one-dimensional PDFs of both E and 
7 are centered on the GR values. This happens both for the 
marginalised cases and when fixing one of the parameters 
at the GR value. This shows that the amplitudes of the S10 
COSMOS cosmic shear measurements agree with those of 
the WMAP 7-year CMB measurements. They both show a 
preference for the same value of as- In other words, given 

1 http://lambda.gsfc.nasa.gov/ 



the strong priors of WMAP the amplitude of the cosmic 
shear models is set, with the prefered ones lying close to 
GR. This can also be seen as an indirect test on the am- 
plitude of systematics. Indeed, if there were systematics 
not accounted for in the data, or in the estimation of the 
redshift distribution, they would make the signal less com- 
patible with the WMAP priors and could show up in the 
results as evidence for modified gravity. Fig. [4] shows the 
PDF of E, marginalised over 7, which corresponds to a la 
constraint of E = 1.02 ± 0.38 . 

The finding that these post-GR parameters are centered 
on the GR values is quite robust, making use of the infor- 
mation on the redshift evolution of the amplitude from the 
tomographic correlations. By repeating the analysis using 
only the wide bin auto-correlations (which contains 60% of 
the galaxies) we found a slight change in the (ctsjE) de- 
generacy, shown in Fig. [5j and a shift on the mean of the 
posterior distribution of E to £ = 0.87 . 

6. Summary and Conclusions 

We used the new COSMOS cosmic shear data to test pos- 
sible deviations from General Relativity. The departure 
from GR is parametrized in a phenomenological model- 
independent way extensively used in the literature. The two 
phenomenological parameters allow to derive two other pa- 
rameters that consistently parametrize the modifications 
to two-point cosmic shear amplitude and growth, and are 
scale-independent and have a scaling redshift dependence 
with the ACDM background. 

A more realistic parametrization should take 
into account three regimes of modified gravity 
(|Hu & : Sawic kT"2007[) . allowing for compatibility with 
the background expansion on large scales and agreement 
with GR on sub-galactic scales. Hu & Sawicki (2007) 
give a fitting formula for the non-linear matter power 
spectrum across the regimes. That formula was con- 
firmed for particular classes of theories by perturbation 
theory f Koyama et al. 2009[ ) and N-body simulations 
( Oyaizu et al. 2008[ ). In particular for f(R), according to 
Koyama et al. (2009) it shows systematic deviations of 
at least 10% from the formula of Smith et al. (2003), for 
scales smaller than k sa 0.8/iM pc _1 . COSMOS has an 
average source redshift of z 1.5 and is mainly sensitive 
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to lensing structures at around half comoving distance to 
the source. Therefore, this limiting scale corresponds to 
6 « 20arcmin and the COSMOS range is affected. We 
do not apply this result to our model-independent study, 
but in general a systematic increase in the amplitudes 
of the models implies that cosmic shear becomes more 
sensitive to E, leading to a more precise determination of 
this parameter for the same data. 

Our results use WMAP priors obtained for ACDM mod- 
els. This is consistent with using a ACDM background 
evolution. The CMB power spectrum is affected by the 
Integrated Sachs- Wolfe effect essentially at the recent epoch 
when the modified gravity effects become important and 
the potentials change significantly over time. This affects 
only the scales large enough compared with the photons 
crossing time. Since the effect is at low redshift those large 
physical scales correspond to very large angular scales. In 
this way, the peak structure is not changed and the pa- 
rameters of modified gravity are not degenerate with the 
standard parameters at perturbation level, allowing to use 
ACDM priors. 

We did not find evidence for inconsistency with GR, in 
agreement with other cosmological tests in the linear regime 
(see |Jain fc Khoury 2010| for a review). In particular, we do 
not find the evidence for deviation reported by Daniel & 
Linder (2010) when using small-scale CFHTLS data. Our 
result using the new tomographic correlations of COSMOS 
corroborates the results obtained by Daniel & Linder (2010) 
using a different COSMOS catalogue ( |Massey et al. 20070 . 
We cannot however make a direct comparison of the con- 
straints, since we assume a different redshift-dependence 
and different priors. 

Our analysis relied on the tomographic amplitude of 
2-pt cosmic shear, and did not allow to break the E — 7 de- 
generacy. The scale-dependence of the modifications, which 
is more distinctive on scales larger than those probed 
by COSMOS, should be explored with future wide high- 
precision weak lensing surveys. Such an analysis will have 
a better chance of reaching our original goal, i.e., to inves- 
tigate if given strong cosmological priors there would be 
room to compensate a modification of the predicted ampli- 
tude of the cosmic shear signal with a modification of the 
growth of matter perturbations. 
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